Search for: All records

We have imaged the entirety of eight (plus one partial) Milky Way (MW)–like satellite systems, a total of 42 (45) satellites, from the Satellites Around Galactic Analogs II catalog in both Hαand Hiwith the Canada–France–Hawaii Telescope and the Jansky Very Large Array. In these eight systems we have identified four cases where a satellite appears to be currently undergoing ram pressure stripping (RPS) as its Higas collides with the circumgalactic medium (CGM) of its host. We also see a clear suppression of gas fraction (M_HI/M_*) with decreasing (projected) satellite–host separation—to our knowledge, the first time this has been observed in a sample of MW-like systems. Comparisons to the Auriga, A Project Of Simulating The Local Environment, and TNG50 cosmological zoom-in simulations show consistent global behavior, but they systematically underpredict gas fractions across all satellites by roughly 0.5 dex. Using a simplistic RPS model, we estimate the average peak CGM density that satellites in these systems have encountered to be$\log {\rho }_{\mathrm{cgm}}/\mathrm{g}{\mathrm{cm}}^{-3}\approx -27.3$. Furthermore, we see tentative evidence that these satellites are following a specific star formation rate to gas fraction relation that is distinct from field galaxies. Finally, we detect one new gas-rich satellite in the UGC 903 system with an optical size and surface brightness meeting the standard criteria to be considered an ultra-diffuse galaxy.

 

The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the circumgalactic medium (CGM) is expected to result in a bow shock that leads the LMC’s gaseous disk. In this letter, we use hydrodynamic simulations of the LMC’s recent infall to predict the extent of this shock and its effect on the Milky Way’s (MW) CGM. The simulations clearly predict the existence of an asymmetric shock with a present-day standoff radius of ∼6.7 kpc and a transverse diameter of ∼30 kpc. Over the past 500 Myr, ∼8% of the MW’s CGM in the southern hemisphere should have interacted with the shock front. This interaction may have had the effect of smoothing over inhomogeneities and increasing mixing in the MW CGM. We find observational evidence of the existence of the bow shock in recent Hαmaps of the LMC, providing a potential explanation for the envelope of ionized gas surrounding the LMC. Furthermore, the interaction of the bow shock with the MW CGM may also explain the observations of ionized gas surrounding the Magellanic Stream. Using recent orbital histories of MW satellites, we find that many satellites have likely interacted with the LMC shock. Additionally, the dwarf galaxy Ret2 is currently sitting inside the shock, which may impact the interpretation of the reported gamma-ray excess in Ret2. This work highlights how bow shocks associated with infalling satellites are an underexplored yet potentially very important dynamical mixing process in the circumgalactic and intracluster media.

 

Low-mass galaxy pair fractions are understudied, and it is unclear whether low-mass pair fractions evolve in the same way as more massive systems over cosmic time. In the era of JWST, Roman, and Rubin, selecting galaxy pairs in a self-consistent way will be critical to connect observed pair fractions to cosmological merger rates across all mass scales and redshifts. Utilizing the Illustris TNG100 simulation, we create a sample of physically associated low-mass (10⁸<M_*< 5 × 10⁹M_⊙) and high-mass (5 × 10⁹<M_*< 10¹¹M_⊙) pairs betweenz= 0 and 4.2. The low-mass pair fraction increases fromz= 0 to 2.5, while the high-mass pair fraction peaks atz= 0 and is constant or slightly decreasing atz> 1. Atz= 0 the low-mass major (1:4 mass ratio) pair fraction is 4× lower than high-mass pairs, consistent with findings for cosmological merger rates. We show that separation limits that vary with the mass and redshift of the system, such as scaling by the virial radius of the host halo (r_sep< 1R_vir), are critical for recovering pair fraction differences between low-mass and high-mass systems. Alternatively, static physical separation limits applied equivalently to all galaxy pairs do not recover the differences between low- and high-mass pair fractions, even up to separations of 300 kpc. Finally, we place isolated mass analogs of Local Group galaxy pairs, i.e., Milky Way (MW)–M31, MW–LMC, LMC–SMC, in a cosmological context, showing that isolated analogs of LMC–SMC-mass pairs and low-separation (<50 kpc) MW–LMC-mass pairs are 2–3× more common atz≳ 2–3.

 

The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW’s stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, versus fuzzy DM (FDM) withm_a= 10⁻²³eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake’s density by ∼10%, and holds the wake together over larger distances (∼50 kpc) than if self-gravity is ignored. The DM wake’s mass is comparable to the LMC’s infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is ∼20% dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars’ kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM versus FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70–100 kpc.

 

The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the “Timing Argument,” which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history and including the recent pericentric passage of the Large Magellanic Cloud (LMC), and recent work has found that the MW disk is moving with a lower bound “travel velocity” of ∼32 km s⁻¹with respect to the outer stellar halo. Previous Timing Argument measurements have attempted to account for this nonequilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10%–12% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume.

 

Abstract We present analysis of the proper-motion (PM) field of the red clump stars in the Large Magellanic Cloud (LMC) disk using the Gaia Early Data Release 3 catalog. Using a kinematic model based on old stars with 3D velocity measurements, we construct the residual PM field by subtracting the center-of-mass motion and internal rotation motion components. The residual PM field reveals asymmetric patterns, including larger residual PMs in the southern disk. Comparisons of the observed residual PM field with those of five numerical simulations of an LMC analog that is subject to the tidal fields of the Milky Way and the Small Magellanic Cloud (SMC) show that the present-day LMC is not in dynamical equilibrium. We find that both the observed level of disk heating (PM residual rms of 0.057 ± 0.002 mas yr −1 ) and kinematic asymmetry are not reproduced by Milky Way tides or if the SMC impact parameter is larger than the size of the LMC disk. This measured level of disk heating provides a novel and important method to validate numerical simulations of the LMC–SMC interaction history. Our results alone put constraints on an impact parameter ≲10 kpc and impact timing <250 Myr. When adopting the impact timing constraint of ∼140–160 Myr ago from previous studies, our results suggest that the most recent SMC encounter must have occurred with an impact parameter of ∼5 kpc. We also find consistent radial trends in the kinematically and geometrically derived disk inclination and line-of-node position angles, indicating a common origin. 

Abstract A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N -body simulations of a spherical MW halo that recently captured a massive (1.8 × 10 11 M ⊙ ) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter particles from the present-day random distribution. Instead, the orbital poles of particles beyond R ≈ 50 kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor δ ρ max = 30% (50%) with respect to underdense regions and δ ρ iso = 15% (30%) relative to the isolated MW simulation (no LMC) between 50 and 150 kpc (150–300 kpc). The clustering appears after the LMC’s pericenter (≈50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: (1) the LMC shifts the velocity and position of the central density of the MW’s halo and disk; (2) the dark matter dynamical friction wake and collective response induced by the LMC change the kinematics of particles; (3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least ≈1:10 mass ratio), causing the clustering of orbital poles of halo tracers. 

Abstract We report the first 3D kinematical measurements of 88 stars in the direction of several recently discovered substructures in the southern periphery of the Large Magellanic Cloud (LMC) using a combination of Gaia proper motions and radial velocities from the APOGEE-2 survey. More specifically, we explore stars in assorted APOGEE-2 pointings in a region of the LMC periphery where various overdensities of stars have previously been identified in maps of stars from Gaia and DECam. By using a model of the LMC disk rotation, we find that a sizable fraction of the APOGEE-2 stars have extreme space velocities that are distinct from, and not a simple extension of, the LMC disk. Using N -body hydrodynamical simulations of the past dynamical evolution and interaction of the LMC and Small Magellanic Cloud (SMC), we explore whether the extreme-velocity stars may be accounted for as tidal debris created in the course of that interaction. We conclude that the combination of LMC and SMC debris produced from their interaction is a promising explanation, although we cannot rule out other possible origins, and that these new data should be used to constrain future simulations of the LMC–SMC interaction. We also conclude that many of the stars in the southern periphery of the LMC lie outside of the LMC plane by several kiloparsecs. Given that the metallicity of these stars suggests that they are likely of Magellanic origin, our results suggest that a wider exploration of the past interaction history of the Magellanic Clouds is needed. 

ABSTRACT Carbon enhanced metal poor (CEMP)-no stars, a subset of CEMP stars ($\rm [C/Fe]\ge 0.7$ and $\rm [Fe/H]\lesssim -1$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $M_{\rm vir}\approx 10^8{\, \mathrm{ M}_\odot }$ and $M_{\ast }\approx 10^3-10^4{\, \mathrm{ M}_\odot }$ at z = 0, as well as in the halo of the Milky Way (MW). These CEMP-no stars are local fossils that may reflect the properties of the first (Pop III) and second (Pop II) generation of stars. However, cosmological simulations have struggled to reproduce the observed level of carbon enhancement of the known CEMP-no stars. Here, we present new cosmological hydrodynamic zoom-in simulations of isolated UFDs that achieve a gas mass resolution of $m_{\rm gas}\approx 60{\, \mathrm{ M}_\odot }$. We include enrichment from Pop III faint supernovae (SNe), with ESN = 0.6 × 1051 erg, to understand the origin of CEMP-no stars. We confirm that Pop III and Pop II stars are mainly responsible for the formation of CEMP and C-normal stars, respectively. New to this study, we find that a majority of CEMP-no stars in the observed UFDs and the MW halo can be explained by Pop III SNe with normal explosion energy (ESN = 1.2 × 1051 erg) and Pop II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $\rm [C/Fe]\gtrsim 2$, corresponding to the absolute carbon abundance of $\rm A(C)\gtrsim 6.0$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $\rm [C/Fe]\approx 3-4$, by adopting faint SNe, it is still challenging to reproduce CEMP-no stars with extreme level of carbon abundance of $\rm A(C)\approx 7.0-7.5$, observed both in the MW halo and UFDs.